Effect of point defects on heat capacity of yttria-stabilized zirconia 
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First-principles calculation and anharmonic dynamical theory were used in sequence to explain a 
large excess heat capacity observed in yttria-stabilized zirconia in comparison with the additive rule 
value. It is found that the excessive shape of heat capacity decays gradually with the Y2O3 doping 
when the number of environmentally different O sites falls to its zero value at 33 mol % Y203-Zr02 
due to the Y atoms adjacent. The model and results presented in this work provide a key insight 
into the complex behaviour and characterization of fast-ion conductors. 
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Point lattice defects both native and artificial, such 
as vacancies, self-interstitials, solute and substitute 
atoms, strongly affect the most fundamental proper- 
ties of materials. In yttria-stabilized zirconia (YSZ), 
the material with numerous commercial applications [2( , 
formed by the addition of Y2O3 to ZrC-2, the triva- 
lent dopant Y 3+ substitute for some of the host cations 
and, in order to maintain charge neutrality, one O va- 
cancy (Vo) must be created for each pair of dopant 
cations. Y2O3 is used to stabilize the tetragonal (t) 
phase of (Zr02)ioo-a: (¥203)2; over the composition range 
2< x <9 mol % and the cubic (c) fluorite phase with 
4< x <40 mol %. The presence of relaxed Vo and Y sub- 
stitute atoms make the local atomic environments of YSZ 
rather different from the stoichiometric high-temperature 
t- and c-polymorphs of pure Zr02 whose 8-fold cations 
are distorted and perfectly coordinated, respectively. In 
YSZ, the average cation coordination number, ranged be- 
tween 7 and 8, is reduced gradually with increasing Y2O3. 
If the Vo associates with Zr ions, as the X-ray absorp- 
tion findings Q suggest, it may support a coordination- 
driven ordering model of stabilization. Q Placing Y in 
the next nearest neighbour (NNN) cation positions al- 
lows the coordination of Zr in the nearest neighbours 
(NN) sites to the Vo to be similar to the monoclinic ZrC>2 
arrangement, with Zr located in 7- fold coordination en- 
vironments, while Y remains 8-fold coordinated. 

The electron energy-loss near-edge structure (ELNES) 
of YSZ demonstrates the features of the experimental O 
iVedge shapes which depend on the crystal structures 
and the Y2O3 composition. |{| In modelling the most im- 
portant and widely used properties of YSZ the relaxation 
mechanism is essential. The effects of relaxation were 
treated within the plane-wave, pseudopotential based 
free energy molecular dymanics (FEMD) technique. [|J 
Using the relaxed configurations (RC) of point defects, 
the ELNES calculations were carried out.jjj Since the 
results of calculation show very good agreement with the 
experimental O -ftf-edge signal, it seems that theory re- 
flects the realities of relaxation. There is much yet to 
be learned about the fundamental nature of transport 
processes in YSZ. Thus, our motivation is to make a con- 
tribution to the microscopic modelling and tailoring the 



heat capacity of YSZ. 

A large excess heat capacity (EHC) has been observed 
in YSZ at low and room temperatures for 7.8, 9.7 and 
11.4 mol % Y2O3 compared to the heat capacity calcu- 
lated from those of pure ZrC>2 and Y2O3 by the additiv- 
ity rule: AC P = C^ sz - [(1-x) C* r ° 2 + x C^° 3 }. The 
shape of AC P (T) is very similar to that of the Schottky 
disorder. The same behaviour of heat capacity may yield 
the two-level system (TLS), quite a simple model, used 
to simulate the low-temperature properties of amorphous 
insulators: C^ LS ^k B (A/T) 2 e - A / T /( e - A / T +l) 2 , where 
A is the energy difference between the levels. C p LS 
shows a maximum at temperature T m ~ 0.42 A when the 
filling-rate variation for the occupied TLS states reaches 
its maximum value. At T > T m , all TLS states are oc- 
cupied equally and therefore heat capacity decreases as 

C TLS „ 1/T 2_ 

It appears that for an anharmonic oscillator (AO), 
which dynamics is defined by the double-well poten- 
tial U(x)=— x 2 +ax 3 +bx 4 , (a,b >0) with a maximum at 
xq=0 and two minima at x\ and X2, separated by the 
energy A=U(x 2 ) — U(xi) >0, the temperature depen- 
dence of heat capacity is similar to that of TLS. The 
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partition function Zao—t^t J dxdv exp[ 
calculated assuming | Uq |3>| Ui |. If U is approxi- 
mated near its minima by the second-degree polynomials, 
U\ ~ (x — xi) 2 u> 2 /2 and U% ~ (x— xi) 2 w|/2 + A where 
u>i are the frequencies of small vibrations at Xi then, us- 
ing the general thermodynamic relations, one can obtain 
= k B £ — ^Pkzhll +k B . In the latter, the 
temperature dependence of the first term, AC (T) , which 
defines the difference in heat capacity between anhar- 
monic and harmonic oscillators, is close to C p LM and fi- 
nally, at wi/o;2=l, AC = Cj LM . Hence, the EHC effect 
may appear as the result of strong anharmonic vibration 
modes either localized or delocalized ones. 

The fact that at x >8 mol % Y2O3 AC P decreases 
7] with increasing the dopant concentration may sug- 
gest that the delocalized anharmonic vibrations cause the 
EHC effect. It is known that the main contribution to 
the structural transformations comes from particular vi- 
brational modes. In pure Zr02, a zone-boundary soft 
phonon, , which breaks the c-symmetry of the O sub- 
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lattice, displacing the O atoms toward their positions in 
the i-phase, may be responsible for the c — t phase trans- 
formation. In YSZ, the experimental data Q are not fully 
clarified for this soft mode because of the static disorder 
in the O sublattice. Using the RC of YSZ, the X^~-like 
phonon was calculated [£| by means of the frozen-phonon 
method for each composition of (Zr02)ioo-:c ^203)3 be- 
tween 3 and 10 mol % Y 2 3 . For 10 mol % Y 2 3 , the 
effective potential has a single minimum and, with de- 
creasing Y2O3 content, the potential develops two min- 
ima. The temperature dependence of the phonon fre- 
quency, calculated within the modified pseudoharmonic 
approximation, quantifies accurately the transition tem- 
perature above which the c-phase is stabilized. In the 3 
mol % Y2O3 case, it may start around room temperature. 
The upper boundary of the single i-phase field in the ex- 
perimental phase diagram |lCj shows a similar rapid drop 
of temperature with Y2O3 concentration. However, the 
X7-like soft phonon of YSZ, in another words, the de- 
localized vibrations can not be responsible for the EHC 
effect because of the potential symmetry. 

Regarding anharmonic effects, the atomic vibrations 
localized near the Vo should be considered as well. YSZ 
is a solid electrolyte in which ionic transport takes place 
by anions moving among their positions by the vacancy 
diffusion mechanism. If any two Vo sites, which an O 
atom occupies before and after a diffusive jump, are non- 
equivalent then theory enables the double- well asymmet- 
ric potential profile. In such a study of diffusion, the 
molecular dynamics technique can fruitfully be used. The 
principal output of the dynamic computation is the set 
of atomic trajectories from which all other properties are 
calculated. Molecular dynamics has been used in the 
past by many authors to calculate the potential energy 
dependant on the position of diffusion atom.pj Even if 
one could calculate the induced trajectories through the 
saddle-point the results were strongly affected by the re- 
action coordinates of all other atoms. The benefits and 
limitations of this molecular dynamics approach have 
been discussed and it is generally accepted now that us- 
ing the large unit cells this scheme is very complicated 
for precise calculation whereas the results, obtained in 
such a way, are not always satisfactory. [TTj 

The features of the computed U profile were demon- 
strated using the single- Vq 11-atom Zr 2 Y207 cell, which 
contains two formula units of ZrC>2 and one formula unit 
of Y2O3, that corresponds to a composition of 33 mol % 
Y203.[I3| The RC, shown in the (a) panel of Fig. 1, was 
used in a further static optimization using the displace- 
ments of 0-2 along the z-direction while all other atoms 
are fixed at their relaxed positions. In this case, labeled 
in the panel (b) of Fig. 1 as a "solid pass" trajectory, the 
0-2 position in the z—0 plane corresponds to the sad- 
dle point with the potential barrier Eb of ~0.2 eV/cell. 
If Zr in the z=0 plane were allowed to relax when the 
0-2 is moved, as shown in Fig. 1 as the "soft pass" 
trajectory, then the equlibrium position of 0-2 is much 
closer to the saddle point while Ej, becomes significantly 
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FIG. 1: Model of Zr2Y2C>7 after relaxation is shown in the (a) 
panel. The plan views shown in (a) represent the views along 
the [001] and [010] directions. The 0-2 displacement along 
[001] is marked by an arrow. In the (b) panel, the curve, la- 
belled as solid pass, shows the change in the FEMD energy 
as a single 0-2 atom is moved from its equilibrium position 
along the z axis while all other atoms are fixed at their re- 
laxed positions. The soft-pass curve, shown after shifting to 
a common zero energy, corresponds to the 0-2 displacements 
when its nearest Zr may relax. 



less. In Zr2Y207, where all O sites are equivalent to 
each other, the Zr/Y environment enables the symmetric 
form of U(x) and therefore can not result in the EHC 
effect. The saddle-point neighorhood plays a key role on 
the instability-barrier formation, which is very sensitive 
to the model parameters. Thus, the dynamics theory of 
ion transport is still very much a scmicmpirical science 
and further simplification, based on the realities for the 
activation energy, should be made. 

Atomistic RC of the (96-y)-atom supercells (y=l,2,...) 
allow to model (Zr02)ioo-:r(Y203) a ;. In the 95-atom cell, 
which corresponds to ~3 mol % Y2O3, the RC with both 
Y atoms in the NNN shell was obtained. , 5] Fig. 2 shows 
the RC energies, expressed as an equivalent temperature 
per 12-atom cell, to give an idea of the relative stability 
and temperature required. In the 94-atom cell, the RC 
with two VoS along the (111) fluorite direction results 
in a 6-fold coordinated and six 7-fold coordinated NN 
Zr atoms. This result supports the interpretation of the 
neutron data by Goff et al. 01- Electrostatic consider- 
ations suggest that VqS should repel. From this point of 
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FIG. 2: Relaxed energies of the 95-atom YSZ cell are plotted 
vs positions of the two Y dopant cations on the NN or/and 
NNN cation shells. The diffusion jump of the lower apical NN 
O allows to construct the double-levels potential U shown in 
the diagram as a sketch. 



view, the tri-Vn RC, with all Vn separated by \/5/2 from 
each other (in units of the c-Zr02 lattice constant a), 
is rather reasonable. The repulsion tendencies between 
vacancies consistent with both electrostatics and the Zr- 
coordination lowering model. Since the shortest distance 
of a/2 between the Vn sites is not energetically favorable 
in YSZ at low Y2O3 concentrations we believe that each 
Vn has six NN O atoms each of which can jump into the 
empty site. 

The intrinsic process of such chaotic O moves should 
occur in YSZ at finite temperatures via the vacancy 
mechanism. In particular, when an apical O in 3 mol 
% Y 2 3 moves into V the RC with the Y NN and NNN 
atoms is created. The difference A=210 K between these 
two RC was used to construct the potential profile, as 
shown in the lower diagram of Fig. 2. If Vn moves con- 
sequently into the next O site of this anion tunnel along 
the z axis, the Y NN-NNNN configuration is realized. We 
take Eb — E a ~1 eV using the activation energy E a of 
YSZ, obtained from the Arrhenius relation for the ionic 
dc conductivity a ~ u>od 2 exp(—E a /ltBT), where wrj 
is the attempt frequency. To form the effective potential, 
characterized by A and the distance moved during a 
single jump d=\ x\ — xi |, the three-parameter polyno- 
mial U(x)= —ax 2 + bx 3 + cx A was used. Since E^ ^ A, 
the O vibrations near the U minima may be estimated 

as oj ~ \ 1 32 ji , where m is the O mass. As a result, 

y md ' ' 

lo w 5xl0 14 Hz is close to that found experimentally [To| 
showing very reasonable model approximation. 

In Fig. 3a the calculated EHC shape AC P {T) = Af TLS '■ 
Cp LS is plotted, with its dependence on temperature 
and dopant concentration. Here Af TLS was obtained in 
such a way that each maximum coincides with the ex- 
perimental EHC amount 7] at 7.76, 9.7 and 11.3 mol % 
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FIG. 3: Calculated AC P for the 7.76, 9.7 and ff.35 mol % 
Y203-Zr02 are shown vs T/A in the (a) panel compared 
with the experimental EHC shape of 7.76 mol % Y2O3. In 
the panel (b), the ratio of JV TLS , obtained using AC P , to 
the number of vacancies N y are shown as □ for each Y2O3 
concentration and connected by eye. 



Y 2 3 . The EHC observations of 7.76 mol % Y 2 3 are 
shown in Fig. 3a to illustrate the similarity All theoret- 
ical AC P shapes, with a maximum at 88 K, have been 
aligned with the experimental peak at T m —7b K using 
the T/A dimensionless scale. The (b) panel of Fig. 3 
shows that the Af TLS /Af y ratio (Af y is the number of Vb) 
falls quickly with increase in Y2O3 content and, finally, 
Af TLS /Afy -^0 at x=33 mol % when the U symmetry for 
diffusion atom should appear. At low concentrations of 
Y2O3, not all possible diffusion jumps lead to the asym- 
metric U . For example, in 3 mol % Y2O3, M y =l whereas 
M TLS falls from six to four since two O NN sites to Vn 
are energetically equivalent due to the Y atoms adjacent. 
The number of such environmentally equivalent O sites, 
which not contribute to the EHC effect, increases with 
increasing the dopant concentration. 

The high-probability O hopping via Vq between the 
low-energy equivalent sites, localized close up to each 
other, can be considered as some sort of the vacancy 
block, which results in the anisotropic absorption peak 
to the internal frictions measured in 9.5 mol % Y 2 O3.0| 
This observation is in agreement with our findings con- 
firming earlier suggestions that some vacancies exist in 
bound states as di-Vn along the [111] fluorite direction. 
As the number of Vb-blocks increases with Y 2 C>3 dop- 
ing, the peak magnitude for the localized relaxation may 
not scale simply with the dopant concentration, reaching 
a maximum at some critical concentration x m correlated 
with that of the EHC maximum. At x > x m , the number 
of sites on the O sublattice with similar local environment 
increases significantly. This may lead to the long-range 
transport of O ions via vacancies reducing therefore the 
local hopping process in Vo-blocks. Regarding the con- 
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centration profile of localized relaxation explained, such 
observations at x > x m may be masked by overlap be- 
tween the isotropic and anisotropic absorption peaks due 
to diffusion and Vo-block relaxation, respectively. 

In summary, our study of YSZ being performed us- 
ing the ab initio relaxed configurations and dynamical 
model of O hopping allows to explain the EHC anomaly 
observed with the Y2O3 doping. In YSZ, as it has been 
demonstrated, the anisotropic absorption peak of inter- 
nal frictions and EHC effect have a common vacancy- 



block nature showing great potential for theory to deal 
with the realistic model of doped ceramics. Because of 
the general physics involved we believe that similar ex- 
cessive behaviour of heat capacity might be observed in a 
wide class of ionic conductors where anion hopping goes 
through non-equivalent positions. 
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